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H . Abstract 

We consider two models of biological swarm behavior. In these models, pairs of 
particles interact to adjust their velocities one to each other. In the first process, 
called 'BDG', they join their average velocity up to some noise. In the second 
process, called 'CL', one of the two particles tries to join the other one's velocity. 
This paper establishes the master equations and BBGKY hierarchies of these two 
processes. It investigates the infinite particle limit of the hierarchies at large time- 
scale. It shows that the resulting kinetic hierarchy for the CL process does not 
satisfy propagation of chaos. Numerical simulations indicate that the BDG process 
has similar behavior to the CL process. 
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1 Introduction 

The derivation of kinetic equations from particle models of swarming behavior has recently 
received a great deal of attention. In biological swarm modeling, the most widely used 
models are particle ones (also known as 'Individual-Based Models') [H [Til CSl El l37] . 
However, to investigate the large scale behavior of biological systems such as fish schools 
or insect swarms, kinetic [31 [101 [ISl [231 [31] and hydrodynamic [HI [311 [3S] models have 
proved to be valuable alternatives. The question of showing a rigorous link between the 
particle and kinetic levels is mostly open. In [5], a mean-field limit of the Vicsek particle 
model [37] is performed and leads to a nonlinear Fokker-Planck equation proposed in 
|17] . A similar program has been performed for the Cucker-Smale model [H [10]. In [2], 
Bertin, Droz and Gregoire propose a binary collision mechanism which mimics the Vicsek 
alignment interaction [37] and formally derive a Boltzmann-like kinetic collision operator. 
There has been no rigorous justification of this derivation so far. The present paper is a 
step in this direction. 

In this work, we investigate two examples of particle systems representative of swarm- 
ing behavior, the so-called BDG and CL processes. These two processes mimic the for- 
mation of consensus in biological groups about the direction of motion to follow. They 
are binary processes. In the first process, called 'BDG' (after Bertin, Droz and Gregoire 
[2]), two interacting particles join their average velocity up to some noise. In the second 
process, called 'CL' (for 'Choose the Leader'), one of the two particles tries to join the 
other one's velocity up to some noise. In this paper, we focus on space-homogeneous 
problems and ignore the spatial variables. Consequently, interactions may happen among 
any pair of individuals in the pool with a certain probability. We also assume that the 
individuals move in a two-dimensional space with unit speed. The state of each particle 
is described by its velocity vector v on the one-dimensional sphere S^. 

The state of an A-particle system can be described by its A-particle probability F/y. 
In the present framework, F/v is a function of the N velocity coordinates (fi, . . . ,fAr) on 
the torus = (S^)^ and of time. The particle dynamics translates into a time-evolution 
equation for Fn called the 'master equation'. In a previous work [8], we have investigated 
the class of 'pair-interaction driven' master equations, of which the BDG and CL master 
equations are members. We have shown that, as A — >■ oo, propagation of chaos holds. A 
propagation of chaos result states that the solution Fjq{t) can be approximated (in a sense 
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to be defined below) by an A^-fold tensor product of tlie single-particle distribution Fi{t) 
provided that this property is true initially. This means that the particles become nearly 
independent and that the system can be described by its single-particle distribution Fi (t) 
instead of the A^-particle distribution F/y. The dimension of the problem is therefore 
considerably reduced. 

To investigate the large limit, it is difficult to work with Fn alone. Indeed, the 
limit of Fn as — ?■ oo is literally a function of an infinite number of variables. The 
functional treatment is simplified by considering the /c-particle marginal -F/v,fc, which is 
the joint probability of any subset of k particles. The number of variables involved in 
i^jv,fc is k and stays fixed as iV — )■ oo. The drawback of this method is that the equation 
satisfied by F^^k depends on the other marginals F^^k' in general. Thus, the equations 
for the (-F/v,fc)fce{i,...,Ar} are all coupled together, forming the so-called BBGKY hierarchy 
[12]. When N = oo, the hierarchy involves an infinite number of coupled equations and 
is called the kinetic hierarchy. Showing a propagation result in the limit N ^ oo involves 
breaking the coupling between the equations in the kinetic hierarchy in some way. 

Consensus formation in swarm models should be associated with the build-up of cor- 
relations between the particles over time. The fact that the BDG and CL models, as 
a result of [8], satisfy a propagation of chaos result is counter-intuitive. The resolution 
of this paradox lies in the investigation of time scales. Indeed, the result of |8] is only 
valid on finite time intervals at the kinetic scale. On this time scale, the number of col- 
lisions undergone by each particle is bounded independently of A^. The present paper 
investigates whether correlation build-up happens at larger time scales. 

Large time scales are investigated thanks to an appropriate time rescaling, i.e. a 
change in the time time. However, the dynamics must also be rescaled in some way to keep 
the leading order terms in the hierarchy finite. Here, the appropriate rescaling consists 
in letting the variance of the noise involved in the process tend to zero. In the BDG 
dynamics, this is not enough. For this reason, we introduce a 'Biased BDG' dynamics, 
where the collision probability depends on the relative velocities of the particles. Then, 
the rescaling also involves a grazing collision limit, i.e. having the collision occur only if 
the relative velocities of the two particles are small. 

The main objective of this paper is two establish the kinetic hierarchies for the rescaled 
BDG and CL processes. We then investigate whether these hierarchies possess solutions 
which satisfy propagation of chaos. For the CL hierarchy, we show that it is never the 
case. In [8], it was already established that the invariant densities (i.e. the stationary 
solutions) do not satisfy propagation of chaos. This was done by looking at the single 
and two-particle marginals only. Here, we extend [S] by showing that the time-dependent 
solution of the CL hierarchy never satisfies propagation of chaos either. We also provide 
a general formula for the /c-particle marginal invariant density. 

Concerning the BDG dynamics, the situation is unclear, in spite of the apparent 
simplicity of the hierarchy equations. We notice that uniform densities are stationary so- 
lutions of the BDG hierarchy. However, the question of uniqueness of stationary solutions 
for this hierarchy is open. There might exist other solutions which do not satisfy the chaos 
property. In a companion paper [9], we investigate the kinetic equation associated to the 
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BDG dynamics. From [8], we know that propagation of chaos is true and that the kinetic 
equation is vahd on the kinetic time scale. The uniform distribution is clearly a stationary 
solution of this kinetic equation. However, in [9], we show that this equilibrium is linearly 
unstable if the noise level is small enough. This suggests the existence of a second class 
of anisotropic equilibria (similar to the Von-Mises equilibria of EI])- The existence 
of multiple equilibria may be a clue that propagation of chaos is broken at infinite time. 
However, these considerations are pure conjectures at this level. 

To improve our understanding, we use numerical experiments. We generate the sta- 
tionary one and two particle marginals by running a large number of independent time- 
dependent runs of the particle dynamics. The experimental results concerning the CL 
dynamics consolidate the theoretical findings. In particular, the theoretical and numeri- 
cal 2-particle correlations show remarkably good agreement. The experimental study of 
the BDG dynamics shows a similar behavior to the CL dynamics. For this reason, it 
should be expected that the BDG dynamics lacks chaos property on the large time scale. 
However, a rigorous result in this direction is not available yet. 

In the literature, propagation of chaos has been mainly investigated in the context 
of the Boltzmann equation and its caricature proposed by Kac. Early works involve 
the names of Kac, Lanford, McKean and others |26l |28l |30]. They have initiated a 
considerable activity [221 [2H |35]. A new approach has been recently developed in |33j . 
Rates of convergence towards chaos have been investigated mainly in the context of the 
Kac model recently in El EDI [23 EH] . 

Lachowicz [27| has recently considered a class of biologically motivated Markov jump 
processes, and proves a propagation of chaos result as part of the derivation of kinetic 
and macroscopic equations. Some of the basic ideas going into his proof, as well as 
into ours, are present in the original works of Kac [SB], but the models, and hence the 
implementation of these ideas, are quite different. 

The outline of the paper is as follows. In section [2l we present the two particle 
processes that we are interested in. In section [3l we derive the master equations and 
BBGKY hierarchies of these processes. Section HI is the core of the paper. It performs the 
iV — )■ oo limit in the rescaled hierarchies and develop the consequences that result from 
it. Section El reports on the numerical experiments. A conclusion is drawn in sectional 
Finally, two appendices collect the technical proofs of the main results of the paper. 

2 Particle models on the circle 

We consider systems of N particles on the circle S^. Each configuration of the system 
corresponds to a A^-tuple uj = {vi, . . . , Vn) G with being the A^-dimensional torus 

= (T^)^. This state space can be seen as corresponding to the velocities of a system 
of mutually interacting swarming agents (see [8] for details). 

The dynamics is first defined as a time-discrete dynamics. Let denote the value at 
the n-th iterate. We now consider different rules for passing from to w""'"^. 
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2.1 The BDG dynamics 



This dynamics is named after Bertin, Droz and Gregoire [2] who introduced it as a model 
for the Vicsek dynamics Given the state a;" G at time n, finding u^~^^ consists 
of the following steps: 



pick an un-ordered pair {i,j) {i < j) randomly according to a uniform distribution, 
i.e. with probability 

Pii = ^^ (2.1) 

''^ N{N-1) ^ ' 



and compute an 'average velocity' 



i,n = A^i^. (2.2) 



Then define 

yn+l=^yn^ ^n+l ^ ^n^n^ (2.3) 

where and are two independent identically distributed random variables on 
distributed according to the probability g{w). The notations use the multiplicative 
group structure of S^. We suppose that g is symmetric: 



g[w ) = g[w) 



2.2 The biased BDG dynamics 

In the sequel, we will consider a 'biased' version of the BDG dynamics defined as follows: 
Let H{w) be a function w ^ E>\ ^ H{w) G [0, 1] be given, where §^ = {2; G | Re 2 > 0}. 
We also assume that H is symmetric: H{w) = H{w*). The biased BDG dynamics is 
similar to the BDG dynamics except for an acceptance-rejection procedure based on H. 
Namely, the procedure consists of the following steps: 

- pick an un-ordered pair (i, j) {i < j) randomly according to a uniform distribution, 
i.e. with probability (12. ip and compute according to (12. 2p . 

- With probability H{v^*v^), perform the collision according to (12. 3p . where again, 

and tf" are independent identically distributed random variables distributed 
according to the probability g{w). 

- With probability 1 — H{y'^*v^), ignore the collision, i.e. define 
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The phase of the quantity v^*v^ is the angle between v^j and v^. The biased BDG 
dynamics performs the colhsion with a probabihty depending on this angle. For instance, 
one can imagine that the collision will be performed more frequently if this angle is small 
than if it is large. It is a straightforward matter to notice that Re(i)^*w") > 0. This point 
will be proved below. This is why the function H needs only be defined on It is 
also easy to see that v^*v'j = {v^*v^)*. Since H is symmetric, the probability H[v^*v^) 
is invariant under exchange of i and j. We also assume that the time unit is chosen such 
that 

maxiJ(u;) = 1. 

In the sequel, BDG will refer to the biased-BDG dynamics except explicitly mentioned. 
2.3 The CL dynamics 

This dynamics is named after the acronym 'Choose the Leader'. It consists of the following 
steps 

- pick an ordered pair (i, j), {i 7^ j) randomly according to a uniform distribution, 
i.e. with probability 

" iV(iV-l)- ^^-^^ 



Define 

where w'"" is a random variables on distributed according to the probability g{w) 



3 Master equations and kinetic hierarchy 
3.1 Master equation 

In this section, we recall the definition of the master equation of the particle system. 
We first introduce the Markov transition operator Q^. It acts on functions $(0;), with 
u G T^, as follows: 

Qn^cu) = E($(a;"+^) | cj" = u), 

where E is the expectation over all stochastic processes involved in the dynamics sending 

oj'' to 

Let now FJ^{u) be the A^-particle probability distribution function at iterate n. Then, 
by definition, FJ^ is such that 

E($(u;")) = I F^{uj)^{u)du. 
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Therefore, by the standard properties of conditional expectations, we have: 



E($(w"+^)) = j F'^+\u)^{u)duj 
= E(E($(w"+i) Iw")) 



= j ^{uj){Q*^F^){uj)du, (3.1) 
where Q*p^ is the adjoint operator to Q^. The dynamics of is thus: 

pn+l r)* rpn 

^ N — ^N-^N- 

To pass to a continuous-in-time dynamics, we assume that the colhsion times for a 
given particle occur according to a Poisson stream with rate z/. Since there are particles, 
the total collision frequency is of the order of Nu. Then, the time-dependent A^-particle 
distribution function F^ = Fjv(ci;,t) satisfies the following master equation: 

^Fn = uN L^Fn, L*^ = Q*^- Id. (3.2) 

We refer e.g. to [8] for details. The weak form of the master equation is given for any 
test function $ e C°{T^) by: 

d_ 

dt Jf 



FN{uj,t)^{u)du = uN / FN{u,t)LN^{uj)du, L^v = Qiv - Id. (3.3) 

AT JjN 



Because the particles are identical, and they cannot be ordered in a natural way, it 
is natural to assume that the initial distribution is invariant under permuation of the 
indices, and this invariance is then preserved by the dynamics. 

The derivations of the master equations for the BDG and CL dynamics are performed 
below. Before this, we recall the notion of hierarchy and propagation of chaos in the next 
section. 

3.2 Kinetic hierarchy and propagation of chaos 

We first define the /c-particle marginal F/v,fc of F^. For k G {1, . . . , A^}, F^^kivi, . . . ,Vk) 
is defined by 



FnA'^i, ...,Vk)= Fn 



(wi, . . . , Ufc, Vk+1, ...,vn) dvk+1, . . . , dvN- 



By permutation invariance, without loss of generality, we can choose to integrate out 
the last N — k variables only and the resulting F/v^ is also permutationally invariant. 
The equation for F/v,a; is found by integrating (13.21) over (w^+i, . . . ,vn). In general, the 
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right-hand side of the resulting equation involves higher order marginals, such as -F/v,fc+i- 
Therefore, the equations for the marginals are all coupled to each other, forming the so- 
called BBGKY hierarchy (see [12] for instance). The BBGKY hierarchy is a key ingredient 
in the investigation of the N ^ oo limit because the number of variables involved in a 
given marginal is fixed. By contrast, the number of variables involved in Fjsf equals and 
goes to infinity with A^, which makes the functional treatment more complex. As long 
as A^ is finite, the BBGKY hierarchy is finite (the number of equations in the hierarchy 
is equal to A^) and does not carry more information than the master equation itself. 
However, in the limit A^ — )■ oo, the hierarchy becomes infinite and is called the Boltzmann 
hierarchy or kinetic hierarchy. There is no master equation any more and the kinetic 
hierarchy is the only object that carries the information about the process. 

Of course, it is desirable to break the hierarchy into a finite number of equations. For 
this purposes, one says that the sequence F/v,a; satisfies the chaos property (or is /-chaotic) 
if there exists a function f{v) of the single variable v such that 

k 

FAr,fc(fi, . . . ,ffc) -> ]^/(wj), as N ^ oc, (3.4) 

i=i 

in the weak star topology of measures, for each G N. This expresses that the A^-particle 
probability F/v approaches a product probability as A^ becomes large, and translates the 
fact that the particles become nearly independent in this limit. 

For a solution of the master equation 03.21) . one can only expect this property to be 
true if at least the initial condition satisfies it. By "propagation of chaos" we mean that if 
(13. 4p holds for the initial data, i.e. that (-F/v)|f=o is /o chaotic for some function /o, then 
for all t > 0, there is a function f{-,t) such that FN{-,t) is /(■, t)-chaotic. In general, the 
rate of convergence in (13. 4p depends on k and on t, and only in particular cases can one 
hope for uniform in time estimates (see [32],[53]). For a similar class of Markov processes 
with applications to biology, Lachowicz has proven L ^-convergence with bounds of the 
form CN~^^'^\ with a function rjiT) could be decaying exponentially fast with T |27j . 

If propagation of chaos holds, then, as A^ — )■ oo, one can replace FN^k{t) by the product 
Y[j=i /(^i! t) the hierarchy and get a closed equation for f{v, t). The resulting equation 
for f{v,t) is a kinetic equation. Combining the BBGKY hierarchy and a propagation of 
chaos result is one of the ways one can derive kinetic equations from N-particle systems 
(see e.g. [28] in the case of the Boltzmann equation or [26] for Kac's equation). 

In the next sections, we derive the hierarchy for both the BDG and CL processes. 

3.3 Master equation and hierarchy for the BDG dynamics 

The following proposition establishes the master equations for the BDG dynamics. The 
master equation for the BDG dynamics (i.e. when H is identically equal to 1) as been 
previously established in [8] . We assume that the Lebesgue measure dv on is normalized 
so that dv = 1. 
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Proposition 3.1 The master equation for the BDG dynamics is given by 



OF. 



N 



dt 



(fl, ...,VN,t) 



2u 



N 



t<j 



2 / H{z) g{u*Vi) g{u*Vj) Fn{vi, . . . ,uz, . . . ,uz*, . . . ,V]\[) dudz 

'{ti,2)G§lx§], 



H{v*.Vi) Fn{vi,...,vn)\, (3.5) 



where uz and uz* are on the i-th and j-th positions respectively. 



The proof of this Lemma is given in Appendix A, section mi We now turn to the BBGKY 
hierarchy and state the: 



Proposition 3.2 Let A; G {1, . . . , A^}. The k-particle marginal F^^k of the solution of the 
BDG master equation \3. 5\] satisfies: 



OF, 



N,k 



dt 



(fl, ...,Vk,t) 



2iy 



N-1 



< 2 / H{z) g{u*Vi) g{u*Vj) FN,k{vi, . . . , mz, . . . , uz* , ...,Vk)du dz 

i<j<k [ Jiu,z)&^xnX 



- H{v*,Vi) FN,k{vi, . . . ,Vk) 



(3.6) 



H{z) g{u*Vi) F7v,fe+i(wi, . . . ,uz, . . . ,Vk, uz*) du dz 
- j ^ H{v*^.^^Vi) FN,k+i{vi, . . .,Vk+i) dvk+i^ . (3.7) 



As examples, we write the first two elements of the hierarchy. For the one-particle marginal 
equation, the first sum is empty and the only remaining term corresponds to the choice 
j = 1 in the second sum. Therefore, the equation is written: 



dF^ 
dt 



{vi,t) = 2u<2 



H{z) g{u*vi) Fn^2{uz, uz*) du dz 
j H{vl2'"i) FN,2{yi,V2) dv^ . 



(3.1 
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For the two-particle marginal, there is only one term from the first sum, corresponding 
to = (1,2) and two terms from the second sum corresponding to the choices i = 1 
or i = 2. This leads to 



-{Vi,V2,t) 



dt 



N-1 



2 / H{z) g{u*vi) g{u*V2) F]\f 2{uz,uz*) du dz 

'(«,z)G§lx§l_ 



+ {N-2)<^2 / H{z)g{u*vi)FN^3{uz,V2,uz*)dudz 

'(u,z)G§ix§i 



{u,z)es^xs] 



H{vlr^Vi) Fn,3{vi, V2, V3) dV3. 

H{z) g{u*V2) F/v 3(^1, uz, uz*) du dz 

- / H{V*2^V2) FN,3{Vi,V2,V3)dV; 

J§1 



As anticipated, the hierarchy is not closed. Each level k requires the knowledge of the 
next level k + 1. If propagation of chaos holds, i.e. if 



Fn,2{Vi,V2) ^ Fn, liVi) Fn, l{v2), 

then, (13.91) can be substituted into (13. 8p and leads to 
dF^A 



(3.9) 



yVi,t) = 2ul2 I H{z) g{u*vi) Fn^i{uz) Fn^i(uz*) dudz 



J^^H{vl2Vi)FN,Mdv2^ FnA^i)^ ■ (3-10) 



This is the kinetic equation proposed in [2]. The question to be investigated is whether 
the approximation (13. 9p can be used. 

It can be seen from (13. 5p that the master equation can be put in the form 



dF, 



N 



2u 



dt N 



i<j 



(3.11) 



where / is the identity and Q*ij) is the following binary collision operator: 

Qli,j)FN{Vi,...,VN) = 

H{z) g{u*Vi) g{u*Vj) F/v(t>i, . . . , uz, . . . , uz* , . . . , V]^) du dz + 



(M,2)e§ix§i 
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Its adjoint given by 



+ {l-H{v*^Vi))<l>{vi,...,VN), 



is a Markovian operator operating on $ through Vi and Vj alone. In [8], the general 
framework of master equations of the type fl3.11l) . called pair-interaction driven master 
equations, is investigated. It is proved that propagation of chaos holds on any finite time 
interval [0,T]. However, it is not known if propagation of chaos holds uniformly in time 
and in particular, if the invariant measure (i.e. the equilibrium -Fjv^oo corresponding to 
= in (13. lip ) is chaotic. The case H = 1 corresponds to the unbiased BDG model 
and has been investigated in [S]. 

Remark 3.1 For computational purposes, the weak forms of the master equation and 
hierarchy are more convenient. The weak form of the master equation liS. 5\) is : 



d_ 
di 



Fn{vi, ...,VN,t) $(wi, . . .,VN)dvi. . .dv 



2u 



N 



N-1 



i<j 



(ui ,...,djv)gT' 

givlivl) 9{v*jv'j) $(^1, w., , Vn) dv- dv'j 



- $(ui, ...,vn) > Fn{vi, ...,VN)dvi... dv 



N, 



(3.12) 



for any continuous test function ^{vi, . . . ,Vn) on . The weak form of the hierarchy 
p. Tp is written as follows: 

/ FN,k{vi, . . . , Vk, t) $(t;i, . . . ,Vk)dvi. . . dvk = — — - 
at Jjk iV — 1 



{i'i,...,i'fc)eT'^ 

9{v*^<) 9i^:,v'^) "^{vi, t;;, ... , Vk) dv', dv'j 



^{Vi, ...,Vk) > FN,kivi, ...,Vk) 



- $(t^i, ...,Vk) ) FN,k+i{vi, ...,Vk, Vk+i) dv. 



fc+1 



dvi...dvk, (3.13) 



for any continuous test function ^{vi, . . . , Vk) on T^. 
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3.4 Master equation and hierarchy for the CL dynamics 

Before stating the result, we introduce some notations: We write (t>i, . . . , £>i, . . . , vn) for 
{vi, . . . , Wj+i, . . . , Vn), i.e. we mean that Vi is absent from the hst. We also define: 

[FN]iivi,...,Vi,...,VN) = F{vi,...,Vi,...,VN)dvi. (3.14) 

Proposition 3.3 The master equation for the CL dynamics is given by 

-^(Vi, ...,VN,t) = 



2u 

iV 



+ [Fm]i{vi, . . . , . . . , Wat) ) - Fn{vi, . . . ,vn) \ ■ (3.15) 



The proof of this proposition can be found in [8]. We reproduce it in Appendix A, section 
17.21 for the reader's convenience. We now consider the BBGKY hierarchy. We have the: 



Proposition 3.4 Let A; G {1, . . . , A^}. The k-particle marginal F^^k of the solution of the 
CL master equation Ii3.15\) satisfies: 

dFN.k , ,x 

— [Vx-, ...,Vk,t) = 



dt 



2u 



N-1 



Yl {\9{v*Vj) (FN,k-i 

i<j<k ^ ^ 



+ F^h-iivi, ...,V, 



+(iv-*)E5|/ . 



(wi, . . . . . . . . . ,Ufc) + 

, • • • , -I, • • • , ^j, . . . , ffe)^ - FN,k{vi, ■■■,Vk) 

givl^^Vi) FN,kivi, . . . ,Vi, . . . ,Vk, Vk+i) dvk+i 



-FN,k{vi, ...,Vk) 



. (3.16) 



As examples, we write the first two elements of the hierarchy (they have been previously 
established in [8]). For the one-particle marginal equation, we get: 



^^^'^-(t;i,t) = ly [ I 9{v2Vi) Fn^i{v2) dv2 - Fn,i{vi) ] , (3.17) 



dt 
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and for the two-particle marginal, we have: 



— [Vl,V2,t) 



dt ' " ' ' N 
2 



^g{vlv2) {Fn^i{Vi) + Fn^i{v2)) - FjV,2(^'l, ^^2) 



+ {N - 2) ( ^ <j / giv^Vi) Fn,2{v2, fs) dV3 + I g{vlv2) Fn,2{vi,V3) dv3 

-Fn,2{Vi,V2] 



By contrast to the BDG hierarchy, the CL hierarchy is closed at any order. This is 
a very remarkable feature of this model, due to the fact that the pair interaction only 
acts on one of the variables. The CL master equation fl3.15p can be put in the frame of 
pair-interaction driven master equations (13.111) . with 

Qii,j)FN{vi, . . .,Vn) = 

= ^9{ViVj) (^[Fn]-{vi, ...,Vj,...,VN) + [FN]i{vi, . . . , Oj, . . . , 

Indeed, its adjoint 
QI^.)^vi,...,vn) = 

]- ( $(wi, ...,Vi,..., zv^, . . . , i^at) + <l>(wi, . . . , 2;wj, . . . , Wj, . . . , vn) \ g{z) dz, 

is a Markovian operator acting through Vi and Vj alone. Therefore, the result of [8] 
applies and propagation of chaos is true on any finite time interval. However, again, it 
is not known if propagation of chaos is valid uniformly in time or breakdowns at large 
times. 

In [5], thanks to the closed hierarchy, an analytical formula for the marginals of the 
equilibrium density F/v,oo is given. It is shown that, if the noise g is properly rescaled with 
A^, the equilibrium density is not chaotic. In sections H] and |5l we revisit this example 
with a special choice of the noise rescaling and we illustrate the loss of chaos numerically. 
This counter-example is not in contradiction with the previous result of [5] because of the 
noise rescaling on the one hand and of of the large time scales on the other hand. 

Remark 3.2 Again, we give the weak forms of the master equation and hierarchy of the 
CL dynamics, which are useful for computational purposes. The master equation 1^3. 15\} 
is given in weak form: 

^ / Fn{vi, ...,VN,t) $(wi, . . .,VN)dVi. . . dvN = TT-T X] / ( 



zG 



i<j "'{''iv,'!'iv)eT 

^ (^^{vi, ...,Vi,..., ZVi, ...,vn) + ^{vi, ...,zVj,...,Vj,...,VN)] g{z) dz 

- $(vi, ...,vn)\ F^ivi, ...,vn) dvi . ..dvN, (3.18) 
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for any continuous test function $(fi, . . . ,vn) on . The weak form of the hierarchy 
/13.16\) is as follows: 

FN,k{vu ■■■,Vk,t) ^{vi, ...,Vk)dvi... dvk = J'^ ^ f 



^{vi, ...,Vi,..., zvi, . . . , Vk)FN,k-i{vi, ...,Vi,...,Vj,...,Vk) g{z) dzdvi 



+7: I ^{vi, zvj, ...,Vj,..., Vk)FN,k-i{vi, . . . ,Vi, . . . ,Vj, . . . ,Vk) g{z) dz dvj 



...,Vk) FN,k{vi, ...,Vk) dvi dvj 



i<k -^(^I'-'^i 



1 

rk-i 2 



K+i,^)GT2 



' {vi,...,Vi,...,Vk)eT' 

. . . , zvk+i, ...,Vk) FN,k{vi, ...,Vi,...,Vk, Vk+i) g{z) dz dv^+i 

- / ^{vi,...,Vk)FN,k{vi,...,Vk)dv. 
for any continuous test function . . . , Vk) on . 

4 Rescaled hierarchies and the hmit N ^ 00 
4.1 Noise rescaling 

Large time scales are not covered by the propagation of chaos resuh of |8]. The goal of 
this section is to investigate whether propagation of chaos is still valid for the BDG and 
CL dynamics at large time scales or not. To do so, it is necessary to rescale the noise 
distribution g (and, in the case of the BDG dynamics, the bias function H). Indeed, if we 
rescale time to large time scales, we simultaneously need to rescale the collision operators 
in order to keep the leading order terms finite. The appropriate scalings of g and H 
correspond to respectively a small noise intensity and grazing collision asymptotics. 

We suppose that the noise probability distribution g{v) depends on a small parameter 
e and we denote it by The parameter e will be linked to in such a way that e — > 
as iV — >■ oo. Similarly, we assume that H = H^{v) depends on e and we introduce the 
probability distribution 

he{v) = ^^, j H,{v)dv. (4.1) 

We keep the assumption that 

max He{w) = 1, Ve > 0. (4.2) 



(3 
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It will be more convenient to introduce the phases of the velocities, i.e. we will write 

V = with e G M/27rZ. (4.3) 

We note that dv = dO /{2tx). We assume that Qs and he are deduced from probability 
distribution functions 7 and 77 defined on M by the following scaling relations: 

Hypothesis 4.1 We assume that 

9e{e) = V^G[-7r,7r], (4.4) 

/i,(e) = —//(-), V^G [-7r/2,7r/2], (4.5) 

where 

r'' , , dx r^^^'^ , , dx 

Is = / — , Ve= V{X> 



J~7T/e 27r i-7r/(2e) ^ 

and where 7 and r] are probability densities on M (for the measures || and ^ respectively) 
which are even and have finite second order moments and r^: 

= [ 7(x) a;^ — < 00, = f vi^) — < 00. (4.6) 



00 27r j-00 



The limit e in corresponds to a small noise intensity limit. In h^, £t represents the 
typical relative velocity at which collisions may happen. With relative velocities larger 
than er, the particles have very little probability to collide. Therefore, the limit e — ?■ 
in he represents a grazing collision limit, in a way similar to the grazing collision limit of 
the Boltzmann equation [T6l [T9] . The hypothesis 14. II can be weakened but since the main 
purpose of this paper is illustrative, we do not seek the broadest generality. We stress 
that the statements given in the following sections are formal. 



Remark 4.1 With hypothesis (4jJ_), we have Tie = 0{e). Indeed, with the assumption 
fi4.^ ), we have 

Heie) = vi-)- 

maX[_^/(2£),,r/(2e)] V ^ 

Therefore, 

T^e ~ CqE, Co 



maxR77 
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4.2 Rescaled BDG hierarchy and the hmit N ^ oo 

We introduce the following operator: 



Theorem 4.1 We assume that e = En is linked to N in such a way that en as 
N ^ oo. We rescale time in such a way that t' = 2'Hs^E%t. We assume that for any fixed 
k, Fj\!^k converges in the weak star topology of measures as N ^ oo towards a probability 
measure F[k], uniformly on any finite time interval. Additionally, we assume that F/v,fc 
remains bounded in C^((M/(27rZ))'^) uniformly with respect to N, and with respect to time 
on any finite interval. Then, F^^] is a solution of the following infinite hierarchy: 

. . . , e,, t) = z/ {a^ - T^) ^(Afc+iF[fc+i])(ei, ...A.---.ekA). (4.8) 

i<k 



The proof of this result is given in Appendix B, section 18.21 Of course, this theorem is 
formal because it is not known it F/v,fc satisfies the assumptions. 

In the limit — > oo, the hierarchy takes the form of an inductive sequence of heat-like 
equations on the fc-dimensional torus. When a = t, the evolution of Fj^j takes place at a 
longer time scale. To find this evolution, one must expand the master equations to higher 
order terms in Such higher order expansions are beyond the scope of the present 

paper. 

If the chaos assumption 

k 

F[k]{9,,...9k) = l[F[,]{9,), (4.9) 
is true, then F[i] = F[i](6',t) satisfies the nonlinear diffusion equation: 

This diffusion equation is of forward type (and thus, well-posed in the classical sense) if 
and only if a > r. If a < t, the diffusion equation is of backwards type and is only well- 
posed for specific initial conditions. In the case a > t, the noise added after the interaction 
(measured by cr) is larger than the typical distance between colliding particles (measured 
by r). Thus, in average, the particles are further to each other after the collision than 
before it. The dynamics is then of diffusive type. Conversely, if a < r the particles are 
in average closer to each other after the collision than before. This dynamics produces 
concentrations at a rate which depends on the regularity of the initial data. However, it 



16 



is not known if propagation of chaos holds for this model and the validity of fl4.10p (even 
in the case r > cr) is subject to caution. 

We can look at the invariant densities, i.e. the stationary solutions of the hierarchy 
(14. 8p . The uniform density 

F[k],e^ = 1, \/ke W, (4.11) 

is an obvious equilibrium solution. It satisfies the chaos assumption, i.e. is a /c-fold tensor 
product of the uniform single-particle marginal -F[i],cq = 1- It is not easy to see if this is the 
unique stationary solution of the hierarchy. The right-hand side of (14.81) cancels functions 
of the form (f){6i — 6k). Therefore, it is tempting to think that invariant densities similar 
to those of the CL dynamics (14.231) (see next section) exist. However, a quick check shows 
that it is not the case, unless for uniform distributions. If the uniform densities (14. lip are 
the unique equilibria of the hierarchy, this could be a hint that propagation of chaos could 
be true for the BDG hierarchy. This would be in marked contrast with the CL hierarchy 
which is examined in the next section. However, the numerical simulations performed in 
section [5] seem to indicate that the two kinds of dynamics have a quite similar behavior. 
Therefore, we cannot make any conjecture whether propagation of chaos holds for the 
hierarchy (14.80 . 

4.3 Rescaled CL hierarchy and the hmit N ^ oo 

Theorem 4.2 We assume that 6 = 6^ is linked to N by 

- (4.12) 



We also rescale time according to t' = t/N. We assume that, for any fixed k, Fj^jk 
converges in the weak star topology of measures towards a probability measure on [0,27r]'^', 
as iV — )■ oo uniformly on any finite time interval. Then, the limit F[k] satisfies the infinite 
hierarchy: 

^ ^ . . . ,9k,t) = ^ ^(^F[k-i]{Oi, . . . ,9i, . . . ,9j, . . . ,9k)+ 



dt 



i<j<k 



+F[k-i]{9i, . . . ,9i, . . . ,9j, . . . , 9k)j 5{9i — 9j) — 2F[k\{9i, . . . ,9k] 



i<k * 



If k = 1, only the second term (in factor of a"^) remains. 



We notice that the assumptions on F^^k are weaker than in the BDG case and actually 
close to be satisfied. Indeed, for any given time, we can extract a subsequence F^^k which 
converges in the weak star topology of measures. What is lacking is some uniform time 
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estimate which would allow the extraction of a single sequence on a whole time interval 
and the uniform weak convergence on this interval. 

The link between e and expressed by f l4.12p is tighter than in the BDG case. This 
scaling allows to keep the largest possible number of non-zero terms in the hierarchy. Each 
level of the hierarchy involves a damped heat equation on the torus with a delta source 
term involving lower order terms of the hierarchy. Like in the finite case, the hierarchy 
is closed at any order. This remarkable feature allows us to show that this hierarchy does 
not satisfy the chaos property (14. 9p . This is expressed in the following 

Theorem 4.3 The CL hierarchy (JjJ^ does not satisfy the chaos property ^.9^ . 

Proof: The equations for the first and second marginals are respectively: 

^ = ^ ^ (4 14) 

^ = (F„(90 + F,,(0,)m - 0.) - 2F„ H- ^ + ^) . (4.15) 

Now, let us suppose that propagation of chaos holds. This means that for all solutions 
of (gUD, the function F[2](^i, ^2, = F[i](^i, t)F[i](^2, must be a solution of (1415|) . 
However, for such a F[2], we have: 

^-^(^1,^2,^) = F[i](ei,t)^(^2,t)+F[l](^2,t)^(^l,t) 







tf^^^^ 
' del 






92F[2] 




[ del 


^ del 



d'Fn 



(^2,t)+F[i](^2,t)^(^l,t) 



This implies that 

(F[i](^i,t) + F[^{e,,t))6{e, - e^) = 2F^^{e,,t)F[^{e,,t) 

which has F[i]{e,t) = for only solution. This shows that in general, the CL dynamics 
does not satisfy the chaos assumption. ■ 

We can precise what the first and second marginals of the equilibrium are. 
Proposition 4.4 (i) The only stationary solution of I!i4-H ) is the isotropic measure 



FnUOi) = 1- (4.16) 

(a) The only stationary solution of ( [^.15] ) corresponding to the one-particle marginal 
is of the form 

Fi2]Adi,e2) = M{ei-e2), (4.17) 
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where A^(6') is the solution of 
Ai{6) has the expression: 



{9)+Mi9) = S{9). (4.18) 



Mie) = ^-^^^^^^ {e'/' + e^'-'^') , ^G[0,27r], wUh a = ^. (4.19) 



Proof: -F[2],eq satisfies 



2 ^9^-F[2],cq d'^F[2],cq 



J-^^ + -^[2],cq = S{9i — 02 



del 86, 

Taking Fourier series, it is easy to see that the Fourier coefficients F[2],cq{ni,n2) such that 
^1+^2 7^ are identically zero. Therefore, -P'[2],eq is the sum of a Fourier series which 
only involves wave-numbers (ni, n2) such that ni + n2 = 0. It follows that -F[2],eq is of the 
form fl4.17p with Ai{6) satisfying (I4.18p . Formula ( I4.19p is a simple calculation which is 
left to the reader. ■ 



Remark 4.2 We again verify that the stationary solution of the CL hierarchy does not 
satisfy the chaos property. Indeed, we have 

F[2]A^l,e2) = M{e, -62)^1 = F[i],e,(^l)F[i],e,(^2). 

That the equilibrium distribution is not chaotic was already proved in f^. Here, we have 
shown that the whole time- dynamics of the hierarchy does not satisfy the chaos property. 
Again, this is not in contradiction to the propagation of chaos result of f^, which applies 
to the unsealed dynamics at a shorter time-scale. 

Remark 4.3 Because of (JjJ^j can denote the scaled noise probability by g^. In JEj, 

the scaling of g^q was defined in such a way that 

lim (A^ - 2){gN[n) - 1) := G{n) (4.20) 
exists and is finite, where ^N{n) is the n-th Fourier coefficient of g^: 

9N{n)= / e-^^^(e)-, neZ. 



Here, assuming that 7 defined by ( [^.^[ ) decays at infinity fast enough, we find that G{n) 
resulting from ^.20^ is finite and given by 



G{n) = —. (4.21) 
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In JEl, the correlation function Ai is found from its Fourier series: 



which, with yields 



1 — G{n) 



for our choice of the noise distribution qn- It is a simple computation to check that the 
Fourier series of Ili4-19\ ) precisely gives ^.22^ , showing the consistency between the present 



results and those of J^. 

We can say more about the marginals of the invariant density, namely: 
Proposition 4.5 The k-particle marginal of the equilibrium Fyi^^f.^ is of the form 

F[k],eq{dl, ■ ■ ■ ,dk) = 0fc-l(^l — dk, ■ ■ ■ , dk-l — dk) (4.23) 

= </.fc_i(^i - ^„ . . . , 9,^1 - 6,, -e,,...,ek-e,), Vz G {l, . . . , A;}, (4.24) 

where (f)k~i{^i, ■ ■ ■ ,^k-i) is 2'iT-periodic in each coordinate and permutationally invariant. 
Moreover, 



Jk-i 



(ei,---,a-i)= [ 0fc_i(ei,...,a-i,ao^- (4.25) 



Proposition 14.51 shows that the /c-particle marginal is just a function of the relative phases 
of the particles with respectively to any one given particle. Therefore, up to the phase of 
this given particle, the equilibrium statistics is perfectly known. 

Proof. We prove by induction that the Fourier coefficients F[fc],eq('^i, • • • ,nk) of [^j^eq are 
of the form 

F[k],eq{ni, ...,nk) = c(ni, . . . , n^) 5(ni + . . . + n^), (4.26) 

with permutationally invariant coefficients c{ni, . . . , Uk). The notation 6{n) stands for the 
Kronecker symbol 5no- Indeed, if this is the case, we can write 

(ni,...,ns.)eZ'= 

ni+...+nfc=0 

= ^(^1' • • • ' ^'^-i' -(^1 + • • • + nk-i)) e^("i(^i-^'=)+-+-^-i(^^-i-^^)) (4.27) 

= Y c{ni,...,ni_i,-{ni + ... + ni + ... + nk),ni+i,...,nk) 

(ni,...,ni,...nj,) G Z*-! 

gj(ni(6»i-6'i)+...+ni_i(ei_i-0i)+n,+ i(ei+i-ei)+...+nfc(6'fc-0,)) 
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Thanks to (14271) . we can write (lOSj) and thanks to IK28\i . we have l^iMl, with 

0fc-l(^l> • • • 7 ^fe-l) = 

= Yl • • • ' ^^-1' -(^1 + • • • + ^'^-i)) e^("^«^+-+"'=-i«'=-i) (4.29) 

= ^ c(ni, . . . ,ni_i, -(ni + . . . + rij + . . . + . . . ,nfc) 

(ni,...,ni,...nfc) eZ^-l 

gi(rii5+...+ni_igi_i+n,+i5i + ...+nfe5fe_i) gg^ 

We note that, thanks to the permutation symmetry of c, the two definitions fl4.29p and 
( I4.30p are consistent. Additionally, the so-defined (pk-i is 27r-periodic and permutationally 
invariant. Finally, f l4.25p follows from integrating fl4.24p with respect to 6k- 

We now prove fl4.26p . Taking the Fourier series of dHSl) (with ^ = smce we are 
interested in the equilibrium, we find: 

^''''"''"'•■■-"''=f(E.,."?) + M^--i).£ 

[F[fc_i]^eq('^i, . . . ,ni, . . . ,ni + rij, . . . ,nk) + 

+F[fc_i],eq(ni, ...,ni + nj,...,hj,..., Uk)] . (4.31) 

Now, we proceed by induction, starting from k = 2. In this case, from f l4.16p . we have 
F[i],eq{n') = S{n) and we get 



2 

f (n?rni)+2 

which is of the form (14.260 . Next, inserting (I4.26P at level k — 1 into (I4.3ip . we get: 



^f'^''^^^"^' • • • '"^^ = 4(y...ni) + k(k-i) 



[c(ni, ...,hi,...,ni + nj,...,nk) + 

+c{ni, ...,ni + nj,...,nj,..., Uk)] 6{ni + . . . + Uk), 



which is also of the form f l4.26p . Therefore, (I4.26P is proven by induction on k, which ends 
the proof of proposition H75l ■ 



The particularly simple form of the CL hierarchy leads to many analytical formulas. These 
formulas can be used to compare the theoretical predictions to numerical experiments. 
This is performed in the next section. 
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5 Numerical simulations 



5.1 Operation mode 

The discrete CL and BDG dynamics (see section [2]) have been simulated. The experimen- 
tal protocol is as follows. For a given number of particles, M independent simulations 
are performed (in the experiments, M = 1000). Each simulation is run until an equilib- 
rium is reached. The detection of the equilibrium is detailed below. Once an equilibrium 
has been reached, one particle is randomly picked and its corresponding 6 is collected and 
added to an histogram. Similarly, a pair of particles is randomly picked and the corre- 
sponding pair {61,62) is added to a 2-dimensional histogram. Then, a new simulation is 
started and the procedure is continued until the M simulations have been performed. The 
histograms of the M samples of 6 or (^1, 62) gives experimental access to the steady-state 
one and two particle distribution functions of the process. 

To quantify if the equilibrium state is reached, three macroscopic quantities are ob- 
served: 



The average velocity: t>" = J2i=i 



- The average direction: uj'^ = -^r^ 

- The order parameter: a" = \v"'\'^ 

It is assumed that the equilibrium is reached when the relative difference i^^*^ (where 
(fi is one of the above defined macroscopic quantities) at two iterates n and n + k (where 
AC is a constant (generally k = 1.2N)) is smaller than a fixed tolerance e -C 1, i.e. 



< e. 



Therefore, the process is supposed to have reached an equilibrium if during k time steps 
the relative change in the macroscopic quantities is small. However a lower bound for the 
number of iterations A/" is set in order to have all the particles interact at least once: 

^>XN (5.1) 

Typically the value A = 3 has been chosen. 

Finally, we use a periodized Gaussian with standard deviation as noise distribution: 

gNi6) = 27rJ2 ^e^^, 6 G [-7r,7r]. (5.2) 

^7vV27r 

We use 

= (5.3) 
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to be consistent with the rescahng of proposition 14.21 Specifically, we have 



cr = 27r7, (5.4) 

where a is the standard deviation defined at (14.61) . 7 is a measure of the intensity of the 
noise if the length of the unit circle is fixed to one. We will use 7 = 1/2, 7 = 1/20, 
7 = 1/200, meaning that we expect that the width of the noise distribution (which is 
equal to 27) be comparable to the length of the unit circle or equal to 10% or 1% of it 
respectively. 

For the BDG dynamics, we use a similar law (15. 2p for with tn related to 7' by a 
similar relation as (15. 3p and 7' to r defined at (14.61) by a similar relation as (15. 4p . 

5.2 Results 

5.2.1 Unbiased BDG dynamics 

The equilibrium one and two particle distributions for the unbiased BDG dynamics are 
presented on Fig. [T]as a function of the number of particles for 7 = 0.02. The histogram 
of the one-particle distribution (left) clearly shows that the distribution is isotropic. In- 
deed, in spite of a relatively large noise level, we cannot distinguish any structure. The 
2-dimensional histogram of the two-particle distribution (right) shows a high level of cor- 
relation along the diagonal, due to the small value of 7. Clearly, the 2-particle distribution 
is not a product of two copies of the one-particle distribution and the chaos property is 
not satisfied. 

5.2.2 CL dynamics 

The equilibrium one and two particle distributions are presented on Fig. [2] as a function 
of the number of particles. Similar conclusions as for the unbiased BDG dynamics can 
be drawn, except for a slightly larger spread of the 2-particle distribution function about 
the diagonal. We also see that the spread is about the same for all particle numbers, 
which confirms the relevance of the choice of scaling the noise standard deviation like the 
inverse square root of the particle number. The statistics of the 2-particle distribution 
function for = 10^ is noisy because the convergence to the equilibrium is very slow and 
the simulation was stopped before the equilibrium was reached. Yet, correlations appear 
with a spread of the two-particle distribution about the diagonal of the same order of 
magnitude as for smaller number of particles. 

Fig. [3] explores the dependence of the one and two-particle distribution functions as 
a function of 7, for a given particle number = 10^. For a large value of 7 (7 = 1/2), 
almost no structure appears in the two-particle distribution function. When 7 is reduced 
to 7 = 1/20, a clear correlation along the diagonal builds up. When 7 is further reduced 
to 7 = 1/200, the two-particle distribution exhibits a delta- like behavior, with a total 
concentration on the diagonal. 
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5.2.3 BDG dynamics 

On figure |H the two-particle distribution law of tlie BDG dynamics is presented for 
Y = 10. Tliis means tliat on tlie average, tlie particles are 10 times closer after performing 
a collision than before. We are in the case of backwards diffusion, with o"^ — < 0. The 
one-particle distribution is not plotted as it looks very similar to the other cases. For a 
small particle number = 10^, two-particle correlations along the diagonal are formed. 
Around the diagonal, a hollow zone reflects the fact that if particles are close they collide a 
lot and tend to correlate their velocities. On the other hand, if they are too far away, they 
cannot collide and correlation does not appear. This hollow zone seems a good indicator 
to measure the ratio if is not known. For a larger number of particles = 10'^ or 
= 10^, the results are more noisy, but this is due to the fact that we have to stop the 
simulation before an equilibrium is reached, due to very large CPU times. Indeed, most of 
the collisions are rejected in the acceptance-rejection procedure. The acceptance-rejection 
procedure should be biased in order to restore reasonable simulation times. This has not 
been done yet and will be the subject of future works. 

5.3 Comparing the numerical and analytical solutions 

Using f l4.19p the analytical solution for the equilibrium two-particle distribution function 
of the CL dynamics can be visualized and compared to the numerical solution. Fig. [5] 
shows this comparison for three values of 7: 7 = 1/2, 7 = 1/20, 7 = 1/200. Despite the 
inherent fluctuations due to the noise and the relatively limited number of independent 
experiments {M = 1000), the agreement appears quite convincing. 

6 Conclusion 

We have considers two models of biological swarm behavior, the BDG and CL models. 
They describe consensus formation in a collection of moving agents about the direction 
of motion to follow. We have first formed the master equations and BBGKY hierarchies 
of these two processes. Then, we have investigated their large particle number limit at 
the large time scale. To this aim, we had to simultaneously rescale the noise variance 
of the processes and make it small in order to keep the leading order of the dynamics 
finite. We have examined the resulting kinetic hierarchies and shown that, in the case 
of the CL dynamics, propagation of chaos does not hold, neither at the stationary state, 
nor in the time-dynamical process. We could exhibit the general form of the /c-particle 
correlations at the stationary state. The study leaves the related questions for the BDG 
dynamics mostly open. However, numerical simulations indicate that the BDG dynamics 
has a similar behavior as the CL dynamics. 

More theoretical investigations are on the way to understand the behavior of the BDG 
dynamics closely. In the case of the CL model, the knowledge of the general form of the 
correlations opens the way for the derivation of a generalized kinetic model. Such a model 
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would replace the standard kinetic model which does not hold because of the violation of 
the chaos property. It would involve a kinetic-like equation for the correlations. 
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(b) N = 10^ 




(c) N = 10^ 




pi.2 pi 3pl.'2 ?pl pi;? pi 3pl<'2 Spl 



(d) AT = 10^ 

Figure 1: Unbiased BDG dynamics with 7 = 0.02 and four values of the particle number 
(from top to bottom: = 10^, N = 10^ N = 10^ = 10^). One and two particle 
equilibrium distributions (left and right respectively). Left: the horizontal axis is 6* G 
[0,27r], the vertical axis is the distribution of 6. Right: the square is the domain of 
{01,62) G [0,27r]^. The color coding indicates the 2-particle distribution. 
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(LtH 




p|.'2 pi :ipl.'2 2pl 



(a) 7V = 




|>l<'2 pi 3pl.'2 2|>k 



(b) iV = 




pi 



(c) N ^ 10^ 




pk.'Z pi ^pli 2p\ a pVl pi 3pl/3 2pl 



(d) TV = 10^ 

Figure 2: CL dynamics with 7 = 0.02 and four values of the particle number (from top 
to bottom: = 10^ = 10^, TV = 10"^, TV = 10^). One and two particle equilibrium 
distributions (left and right respectively). Left: the horizontal axis is 6' G [0,27r], the 
vertical axis is the distribution of 9. Right: the square is the domain of {9i, 62) G [0, 27r]^. 
The color coding indicates the 2-particle distribution. 
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pV2 pi 3plj? Spi a ptl2 pi SpL'? 2pl 



200 

Figure 3: CL dynamics with = 10'^ and three values of the parameter 7 (from top to 
bottom: 7 = 1/2, 7 = 1/20, 7 = 1/200). One and two particle equilibrium distributions 
(left and right respectively). Left: the horizontal axis is ^ G [0, 27r], the vertical axis is the 
distribution of 6. Right: the square is the domain of {01,62) G [0,27r]^. The color coding 
indicates the 2-particle distribution. 
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(a) N = 10^ 



(b) N = 10^ 




pi/2 pi 3pi'2 2pi 



(c) N = 10-* 

Figure 4: BDG dynamics with 7 = 0.02 and 7' = IO7 and three values of the particle 
number (from top to bottom: = 10^, = 10^, = 10^). The two particle equilibrium 
distribution only is plotted. The square is the domain of (^1,6*2) G [0,27r]^. The color 
coding indicates the 2-particle distribution. 
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(c) 



Figure 5: Perspective view of the equilibrium two-particle distribution of the CL dynamics 
as functions of {01,62) G [0,27r]^, with (a) 7 = | ; (b) 7 = ^ ; (c) 7 = Left: 
numerical solution. Middle: analytical solution. Right: Difference between the numerical 
and analytical solution. 
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7 Appendix A. Establishment of master equations 
and hierarchies 



7.1 BDG master equation and hierarchy 

Proof of Proposition 13.11 The Markov transition operator for the BDG dynamics is 
clearly given by: 



H{v*jVi) / ^{vi,...,wVij,...,zvij,...,VN)g{w)g{z)dwdz 

+ (1 - H{v*jVi)) ...,Vi,...,Vj,..., vn)^ , 

where wvij and zvij are on the i-th and j-th positions respectively. Here, we recall that 
Jgi dv = 1. We then compute: 



/$ := / QN^ivi, ...,vn) F^ivi, . . . ,VN)dvi. . . dvN 

J (di,...,i)jv)gT" 

H{v*-Vi) / . . . , wVij, zVij, ...,vn) g{w) g{z) dwdz 

+ (1 - H{y*-Vi)) . . . , Wj, . . . , . . . , vn)^ Fn{vi, . . . ,vn) dvi . . . dv^ 

The change of variables {w,z) G — )■ {v[,v'j) G defined by v[ = wVij, Vj = zVij is of 
jacobian unity. Therefore, we deduce that 



-T77^-r.J2 \ I [ H{v*^v,)^vi,...,vI,...,v'j,...,vn) 



aiv-jv'i) giVijv'j) FN{vi,...,VN)dvi... dvN dv[ dv' 



+ / (1 - H{v*jVi)) ...,vn) Fn{vi, . . . ,vn) dvi . . . dvN 

J [vi,...,vm)&T" 

Using (13. 2p leads to the weak form fl3.12p . 
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Now exchanging the notations {vi,Vj) and in the first term, we have: 

= WJ^^^l I I H{v-*v'^<^>{Vi,...,Vi,...,Vj,...,VN) 

g{Vi*Vi) g{Vi*Vj) Fn{vi, . . . ,v'i, . . . ,Vj, . . . , vn) dvi... dvN dv[ dv'j 



{vi,...,VN)eTS'" 



(1 - H{v*jVi)) ...,vn) Fn{vi, ...,VN)dvi... dvj^ 



with 



^^ + 

V- — 



With (O), we deduce that 

H{v[*v[) g{Vi*Vi) g{v'*jVj) Fn{vi, . . . , w^, . . . , w^, . . . , vn) dv[ dv^ 



+ (1- H{v*m))Fn{vi,...,vn) 



U = v'iA, z = v*,v' 



Now, we change variables {y[,v'j) G to {u,z) G S"*^ x §^ such that 
where SJ',_ = {2; G | Re^ > 0}. The restriction z G SJ',. comes from the observation that 
Ke{v*jV[) is non-negative. Indeed 

and since = \v'j\ = 1, we have |Re(f^*fj')| < 1 and 1 + Re(f^*w^) > 0. Reciprocally, it 
is easy to see that for (u, f ) G x there exists a unique pair {v[, Vj) G such that 
u = v[j, z = v*^v[. Indeed, clearly, v[ = uz, and we have uz*v'* = 1, which gives v'j = uz*. 
Finally, using the phases, one sees that dv'^ dv'j = 2du dz. Using this change of variables, 
we have: 

H{z) g{u*Vi) g{u*Vj) F^{vi, . . . , uz, . . . , uz*, . . . , V]\r) du dz 

+ {1- H{v*jV,))Fn{vi,...,vn) 

With definition (13 .2^ for the master equation, we are led to (\3.5h . which ends the proof of 
Proposition 13.11 ■ 
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Proof of Proposition 13.21 We start with the weak form (I3.13p . We take $ only 
depending on (t>i, . . . , Vk) in f l3.12p . We get 



d_ 

di 



Fn{vi, ...,VN,t) ...,Vk)dvi... dvN 

fN 



d f 



,Vk)dvi. . .dvk 



= Ti + T2 + %i (7.1) 

where 7i (resp. T2 resp. Ts) collects the terms at the right-hand side of (13.121) corre- 
sponding to i < j < k (resp. i < k < j ; resp. k < i < j). 

We first note that Ts = 0. Indeed, with k < i < j, since $ only depends on (t>i, . . . , Vk) 
and g is a. probability, the curly bracket at the right-hand side of (I3.12p becomes: 

g{v*^v'i) g{v*jv'j) <l>{vi, ...,Vk) dv[ dv'j - <^{vi, ...,Vk)= 0. 

Concerning the first term 7i, we can just integrate {vk+i, . . . ,vn) out of Fn and the 
result is a just the first term at the right-hand side of (I3.13p . Now, we focus on 72. By 
permutation invariance, all terms of the sum over j are identical. We can collect them 
into one single term for say j = k + 1. There are N — k such terms. We can then integrate 
(f • • • , vn) out of F/v leading to an expression involving only Fn^^+i- Now, for the same 
reason as above, we have : 



/ 



aiv-k+i^'i) 9{v*k+iv'j) Hvi, ...,v'i,...,Vk) dv[ dv'f^ 



k+l 



Collecting these observations leads to the second term of at the right-hand side of (I3.13p . 

Now, we turn to the strong form (13. 7p . The first term at the right-hand side of (13. 7p 
is obtained from the corresponding term of (I3.13P in exactly the same way as in the proof 
of proposition 13.11 We skip the details. We now focus on the second term. We exchange 
notations between Vi and v'^ and we change the notation v^+i into We get 

/ H{v*i^^-^Vi)g{v*i^^^v'i) $(ui, . . . , w,', . . . , Wfc) 

FN,k+iivi, ■■■,Vi,...,Vk, Vk+i) dvi dvk+1 dv[ = 
H{y'il+A)g{v\l^iv,) $(wi, . . . , Ui, . . . , Wfc) 

FN,k+i{,vu . . . , . . . , Wfc, Wfc+i) dvi dv'f,^^ dv[, 
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and changing to the variables (f^iffc+i) € to {u, z) G x such that u = "y^fc+i, 
z = v'il+iv'i, we get 

H{v*k+iv'i)g{v*k+iVi) FN,k+i{vi, . . .,Vi,. . . ,Vk, w^+i) dv[^^ dv[ = 



I 



H{z)g{u*Vi) FN,k+i{vi, ...,uz,...,Vk, uz*) du dz. 



Inserting this expression into the second term at the right-hand side of (13.131) leads to the 
corresponding term of fl3.7l) and ends the proof. ■ 



7.2 CL master equation and hierarchy 

Proof of Proposition 13.31 The Markov transition operator for the CL dynamics is 
clearly 

...,zvj,...,Vj,..., vn)) g{z) dz. 

We then compute: 

Qn^{vi, ...,vn)Fm{vi,..., Vn) dvi . . . dv^ 



{vi,...,vr^)eT" 

. . .,zv^, ...,Vj,.. .,vn)) g{z) Fn{vi, ...,VN)dz dvi . ..dv^ 

This leads to expression (13.181) . 

By the change of variables v'^ = zvi in the first term and = zvj in the second one, 
we get 



-y- 



N(N-l) ^ 2 



(i;i,...,{)j,...,?Jiv)GT"~i >/d'gS1 



^{Vi, ...,Vi,...,Vj,...,VN) giv*Vj) 

[FN]j{vi, ...,Vj,..., VN)dvi ...dvj... dvN dv'j 



+ / '^{vi, . . . ,v[, . . . ,Vj, . . . ,vn) g{v*v'i) 

[F^Yiivi, ...,Vi,..., VN)dvi ...dvi... dvN dv[ \ . 
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By exchanging the roles of the primed and unprimed variables, we get: 



2 f \ 

Wat — TtZI / ^{vi,...,VN)-{g{vlv,j)[FN]-{vi,...,Vj,...,VM) + 

Therefore, using the symmetry of we find: 



+ [FN]iivi, ...,Vi,.. .,vn)} 

In view of (13. 2p . we find (13.151) . which ends the proof of Proposition 13.31 



Proof of Proposition 13.41 We start with the weak form (I3.19p . We take $ only 
depending on (f i, . . . , Vk) in (I3.18p . We get a similar expression as (17. ip where 71 (resp. 
72 ; resp. 73) collects the terms at the right-hand side of (I3.18P corresponding to i < j < k 
(resp. i < k < j ; resp. k < i < j). Again, 73 = 0. Indeed, for A; < z < j, since $ only 
depends on (wi, . . . , Vk) and g is a probability, the curly bracket at the right-hand side of 
(I3.18P becomes. 



^{vi, ...,Vk) g{z) dz - $(t;i, . . . , ffc) = 0. 



We can integrate (t'fc+i, • • • , vn) out of Fat in 7i and get: 



2u 



N -I 



(?)i,...,?)fe)eT'= 



E 



i<j<k '-"'^GS 



1 

.1 2 



zvi, ...,Vk) + $(wi, . . . , zvj, ...,Vj,...,Vk)] g{z) dz 



$(W1, ...,Vk) \ FN,k{Vl, ...,Vk) 



dvi . . . dvk, 



However, we can integrate once more the first term in the curly bracket because $(f 1, 
. . . , zvi, . . . , Vk) does not depend on Vj and similarly with i and j exchanged in the second 
term. By permutation invariance, this can be expressed in terms of -F/v,fc_i. This leads to 
the first term at the right-hand side of (I3.19p . 

By permutation invariance. To, is the sum of N — k copies of the term obtained by 
choosing i = k + 1. We can integrate {vk+2, ■ ■ ■ , vn) out of Fn and obtain an expression 
involving F^^k+i only. Now, since i < k < j and $ only depends on (t>i, . . . , Vk), the curly 
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bracket at the right-hand side of (I3.18P becomes, 



/ {/ i 



^{vi, ...,Vi,...,Vk) + . . . , zVk+1, Vk)j g{z) dz 

-^{vi,...,Vk)^ FN,k+i{vi, ■■■,Vk, Vk+i) dvk+1 
^{vi,...,zVk+i,...,Vk)g{z) dz 

- ^{vi, Vk)^ FN,k+i{vi, ...,Vk, Vk+i) dvk+i. (7.2) 
Collecting these observations, we can express T2 as follows: 



r 2^ 



N -I 



(Di,...,i)(.)eT 



{N-k)y2[ \\[ ^{vi,...,zvk+i,...,Vk)g{z)dz 

dvi . . . dvk, 



i<k '""'+- 

- . . . FN,k+i{vi, . . . ,Ufc+i) dvk+i 

However, we can integrate once more the terms in the curly bracket since . . . , -zt^fc+i, 
. . . , Vfc) does not depend on Vi and similarly $(f 1, . . . , w/;) does not depend on Vk+i- The 
resulting expression involves only F^^k and leads to the second term at the right-hand 
side of (13391) . 

To establish the strong form (13.161) from the weak one (13.191) is easy. We perform the 
change of variables Vj = zvi (resp. Vi = zvj ; resp. Vi = zvk) in the first integral of the 
first sum (resp. in the second integral of the first sum ; resp. in the first integral of the 
second sum) and use the sjTiimetry oi g. ■ 



8 Limit — > cxd in the rescaled hierarchies 
8.1 Preliminaries 

In this section, we state the following straightforward lemma, which will be used through- 
out the proofs below. 

Lemma 8.1 The following estimate holds in the sense of distributions: 
ge = 5 + + 0(5^), K = 5 + + 0{e'). 
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More precisely, assuming that the test function $ G C'^(M/(27rZ)), we have: 



(0)) 



< Ce||<l>||3,oo, 



and similarly for with a replaced by t and the integral is taken on the interval [— 7r/2, 7r/2] . 
We denote by C^(]R/(27rZ)) the space of 2tt periodic, three-times continuously differen- 
tiable functions on M and by ||$||3 oo the supremum of all the derivatives of^ up to order?). 



We will also use the following operators: 
A, 



D. 



'' J 

Q2 Q2 



Q2 



_d d_ 



1.2) 



in addition to Dij defined at eq. (14. 7p . 



8.2 Limit — > oo in the rescaled BDG hierarchy 



Proof of Theorem l4.1i We first write the weak form of the hierarchy fl3.13p . introducing 
the phases (14 ■3p and the rescaled noise distribution and biasing function (14. ip . as 
well as the rescaled time (which will still be denoted by t for the sake of simplicity). Now, 
we use the additive group notations of the phases for the multiplication of two elements 
of §^ We have 



d_ 

dt 



N,k\Oi, . . . , 



efc,t)<i>(^i,...,efc; 



dOi... dOk 



V 



i<j<k 



2 ■ 



(e'„9')e[o,2 



dei^i - dij) dei^'j - dij) ^{Oi, . . . ,9'^, . . . ,9'p . . . ,9k) 777-^2 



d9'i d9'j 



$(^1, . . . ,9k) > Fj^^k{9i, ■ ■ ■ ,9k) 



+ {N-k) 



efc+ie[o,27r: 



i<k 



2 



d9' 

gM-0^k+l)m,■■■,0'„...,9k)-^ 

/vr 



$(6'i, . . . ,9k) > F/v^fc+i(6'i, . . . ,9k, 9k+i. 



d9, 



fc+i 



2tt 



d9i... d9k 

{2ti)^ 



1.3) 



for any continuous test function $(^1, . . . , ^fc) on [0, 27i]''. In this formula, 9i^ = {9i — 9j) /2 
or (9i — 9j)/2 + 7i with — 7r/2 < 9^^ < 7r/2 (modulo 27r) is the phase of the quantity v*jVi, 
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while 9ij is the phase of Vij, i.e. 6ij = {6i + 6j)/2 or {6i + 6j)/2 + Ti and 6i — 6ij = 9ij — 9j = 
e [-7r/2,7r/2] (modulo 27r). 

We first investigate the asymptotics of (18. 3p when e — )■ with fixed A^. Take $ G 
C^(M/(27rZ))'^. We first prove the following formula: 



d f 

7[0,27rl'= 



dOi... dOk 
{2-Kf 



V 



N- 1 



+ {N-k){ I V {cr\A,k+i<^) F^^k+i + (I) ife+l-f7V,fc+l) 

(ei,...A-+i)e[0,27r]'=+i 



(2vr; 



fc+i 



•4) 



Moreover, thanks to the assumptions of the theorem, the o(e^) terms are estimated by 
Cf^ll'^'lls.oo with C independent of on any finite time interval T. 

Proof of ( 18. 4h . We start by considering the first sum in (18. Sp . where i < j < k. We use 
Lemma 18.11 and get: 



/ 



)€{R/{2ttZ))' 



. . . , . . . , . . . , ge{9[ - 9s{9', - %) g ^ 



2 [d9f^ 89^ ' 



7]^, . . . , ^jj, . . . , ^jj, . . . , P 



.9k) + o{e'). 



Now, using this formula, we compute the integral: 



hJ9^ 



(6l,,6lj)e(M/(27rZ))2 



r//9' d9' 

$(^1, . . . , . . . , . . . , 4) ^7,(e: - 9,j) ge{9'^ - ^ ^ 



Af.fc 1,(71, . . . , 



(i^j d9j 
^^'^^2^2^- 



The change of variables {9i,9j) e (M/(27rZ))2 ^ {9,j,9i^) e M/(27rZ)x] - n/2,7v/2], is 
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one-to-one onto and we have 



K.. 



,-6M/ 27r2 



,G[-7r/2,7r/2] 



0ijeK/(27rZ) g[-7r/2,7r/2] 



($F/v^fc)(6'i, . . . , 6'ij + . . . , i/jj — . . . , 



QiA — Q i- 



71 Zn 



Using again Lemma [8.11 we find: 



(GM/ 27r 



$(^1, . . . ,9ij, . . . , 9ij, . . . ,9k 



. . . , 9ij , . . . , ^jj , . . . , 9k ; 



271 



/GM/ 27rZ 



,GM/ 27rZ 



q2 ^ ^ d9- 

-QQ^[FN,k{9i, . . . ,9ij + 9, . . . ,9ij - 9, . . . ,9k,t)]g=Q 

-^[{^FN,k){9i, . . . ,9ij + 9, . . . ,9ij - 9, . . . , 9k)]0=o-^ 



+o(e' 



Using definitions (18. ip and (18. 2p . we get: 

^2 



,GK/{27rZ) 



a 



— {Aij^)FM,k- 



We can write 



6l,jGR/(27rZ) J6»,_j G[-7r/2,7r/2] 



1, . . . , 9ij, . . . , 9ij, . . . ,9k,t) + 0(5^) 



2 2 

y(A,,$)F;V,fc + y[$A,F7V,A: 



I, . . . ,9i, . . . ,9j, . . . ,9k,t) 5{9i^) ^ — h o{e^^ 



2 27r TT 

Going back to variables (6'j, ^j), we find, since 5{9i^) = 26{9i — 9j): 



r 
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This gives the first term of (18. 4p . The second term is found using exactly the same type 
of computation. The details are omitted. This ends the proof of (18.41) . 

End of proof of Theorem 14.11 Now, using Green's formula, we get either the weak 
formulation of (14. 8p : 



d_ 

9t -/[o,27r]fe 
V 



...,^fc,t)<l>(^l,...,^, 



{2-kY 



N -I 



(ei,...,0fe)G[O,2^]fe .<^.<^ 

dO df) 



'{ei,...A+i)e[o,2^]fc+i .<fc 

dO df) 



(27r)'=+i 



or the strong form: 



dt 



'-{9i, . . . , 6'fc, t) 



N-1 



I i<k Je,+^e[o,2n] 



[(Afc+l-^Af,A:+l)<^(6'j — Ok+l) — F]y^k+lDik+l{S{9i — 9k+l))] } 



de 



k+l 



271 



,(8.5) 



.(8.6) 



Now, letting — oo and e = En ^ 0, we find that the leading order term in (18. 5p 
or (18. 6p is the second sum because of the factor N — k. Thanks to the hypotheses of the 
theorem, the other terms in the weak form (18. 5p are either o(e^) or 0{1/N) or o{e\j/N) 
multiplied by ||$||3,oo and constants which are independent of and of time. Therefore, 
passing to the limit in the weak form or equivalently in the strong form (18. 6p in the sense 
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of distributions is allowed and leads to 



k+l{Flk+l]S{9,-9k+l)) + 



dt 

i<k 

W [(Afc+iF[fc+i])5(^. - 4+1) - F[fe+i]Afc+i(5(^. - e^+i))] }^. (8.7) 
It remains to show that 

efc+ie[o,27r] 

= {a^ - r^){D,k+iFik+i]){9i, . . . ,9,, . . . ,9^,9,). (8.8) 

We have: 

\k+iiF[k^i]5{9i — ^fe+i)) = (Aifc+iF[fc+i])(5(^j — 6'fc+i) + 2F[k+i]5"{9i — 9k+i) + 

where 5' and 5" denote the first and second derivatives of the Dirac delta. Using Green's 
formula, we deduce that 

d9k+i 



/ ^ik+iiF[k+i]6{9i - 9k+i)) 

+ -7^) + + 2-— — 



-K?. — ) ) F^k+i]{9i, . . . ,9i, . . . ,9k, 9i) 

O^k+l J 



'd9f 391^,' d9l^, d9k+i'd9i 39^+1 
= {F>ik+iF[k+i\){9i, . . . ,9i, . . . ,9k, 9i). (8.9) 

Now, we clearly have 

, d9k+i 



/ {Dik+iF[k+i])5{9i — 9k+i)— — = 

/ 92 92 92 \ 

On the other hand, since Dik+i{5{9i — 9k+i)) = 45" (6*^+1 — 9i), Green's formula leads to 
I F[k+i]Dik+i{S{9i - 6'fc+i))— ^^ii = 4:—^^{9i, . . . ,9i, . . . ,9k,9i). 



But, by permutation invariance, we have 



rj^o v^i7 ■ ■ ■ i9i, . . . ,9k, 9i) — — —2 — (6^1, . . . ,9i, . . . ,9k, 9i). 
991+1 99f 
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Therefore, 

/ {iDik+iF[k+i])5{9i — 9k+i) — F[k+i]Dik+i{5{9i — ^fc+i))} = 

= —Dik+iF[k+i]{Oi, . . . ,9i, . . . ,9k, 9i). (8.10) 

With ([HSD and flHTTOj) . we easily verify ([83]). Now, inserting ^B) into (Q leads to 
and ends the proof of Theorem 14.11 ■ 



8.3 Limit — > oo in the rescaled CL hierarchy 



Proof of Theorem 14.21 We first write the weak form of the hierarchy (13.19p . introducing 
the phases (14. 3 p and the rescaled noise distribution Qf,, as well as the rescaled time (which 
will still be denoted by t for the sake of simplicity). We have 



d_ 

di 



/ FN.k{9i, . . . ,9k,t) ^{9i, . . . ,9k) 



d9i... d9k 

(27r)*^ 



2u 



{N - 1)^2 



E 



.A-)G[0,27r]'= 



^9,,■■■,9i,...,9, + if,...,9k)F^,k^^{9,,...,9,,...,9,,...,9k)9e{^)^^^ 
$(6'i, . . . ,9j + (f, . . . ,9j, . . . , 9k)FN,k-i{9i, . . . ,9i, . . . ,9j, . . . , 9^) Qei'^) J 



$(6*1, . . . ,9k) FN,k{9i, ■ ■ ■ ,9k 



(e„ej)G[0,27r]2 



d9i d9j \d9i...d9i... d9j ...d9k 
"(2^ I (27r)^'-2 



+ {N-k)Y, [ 

i<k JiOi,-A, 



,A)6[0,27r]fe-l 2 



$(^^1, . . . , 9k+i + ^,...,9k) F^49,, ...,§„..., 9k, 9k+,) 9ei^) ^^^^ 



- /" <f(^i,...,^fc)Fj 

^6lie[0,27r] 

Using Lemma [8.11 we have: 



N,k[^i^ 



. . . ,9k 



d9i...d9i... d9k 



(27r)^-i 



(8.11) 



ipe[0,27r]- 



^9i,...,9i,...,9, + ip,...,9k)ge{^)^ 

Zn 



fr^eS ^2^ 
2 



/ 



9i, . . . , ^i, . . . , ^j, . . . , 6^^) + 0(5^) 



d9i 



2 55)2 I v^i, • • • , 61,, . . . , 6'j, . . . , 6^^) 6{9j -^i)^ + 
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and similarly with i and j exchanged. In the same way, we have: 



(/3G[0,27r] 



$(^1, . . . , Ok+i + ip,...,ek) ge{v) TT 

Ztt 



8ie[0,2n] 



2 dof 



9,,...,e,,...,ek)5{e,-ek^,)-^ + oie' 

Zn 



Moreover, the 0(5"^) term is of the form C£:^||$||3^oo where C is independent of A^. There- 
fore, we have: 



d f 



d9i... d9k 



1v 



(iV- 1)£2 



E 

.i<j<k 



(ei,...A-)e[o,27r]fe 



FN,k-i{9i, ■ ■ ■ ,9i, . . . ,9j, . . . ,9k] 



+F'N,k-i{9i, . . . ,9i, . . . ,9j, . . . ,9k) 



89 



j{9i, . . . ,9k) F/v^fc_i(6'i, . . . ,9i, . . . ,9j, . . . , 9k)+ 



3 



d9f 



^1, • • • 7 9k) FN^k~i{9i, . . . ,9i, . . . ,9j, . . . ,9k 



S{9, - 9i 



— ^i9i, . . . ,9k) FN,k{9i, . . . ,9k 



+ iN-k)J2{[ 

^<k 1-^(^1 



d9i... d9k 
(27r)^ 



6lfc+ie[0,27r] 



efe)G[0,27r: 



1 

k 2 



2 d9f 



{9i, . . . , 9k) FN^k{9i, . . ■ ,9i, . . . ,9k, 9k+i) S{9i — ^^.+1 

d9i... d9k 



d9^ 



k+l 



$(6'i, . . . ,9k) Fj^^k{.9i, . . . ,9k 



(27r)^- 



27r 
.(8.12) 



where, since Fj^^k is a probability, the 0{e^) terms are, as previously, of the form Ce^H^Hs^oo 
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and are uniform in and time. However, by permutation invariance, we have 



$(6'i, . . . , 9k) F^^kiOi, ■ ■ ■ ,0i, ■ ■ ■ ,0k, dk+i) S{9i — 9k+i, 



(ei,...A)G[0,27r]'' 



(0l,...A)e[O,27r]fc W0fc+ie[O,27r] 

— $(6'i, . . . ,9k) F^^ki^i, ■ ■ ■ ,9k) 
$(6'i, . . . , 9k) FN^k{9i, . . . ,9i, . . . ,9k, 9i) 

— $(6'i, . . . ,9k) FN^k{9i, ■ ■ ■ ,9k) 
^$(6'i, . . . ,9k) Fjy^ki^i, . . . ,9k) 

— $(6'i, . . . ,9k) F]^^k{di, ■ ■ ■ ,dk) 

Therefore, fl812D yields 

d9i... d9k 



d9k+i 
2ii 



d9i... d9k 

(27r)^ 



d9i... d9k 



,A)6[0,27r]'= 



(2vr) 



d9i... d9k 

{2ti)^ 



d f 

/ FN,k{9i, . . . ,9k,t) ^{9i, . . . ,9k 

7fO,27rF 



(2vr) 



2u 



E 

.i<j<k 



(0l,...A)G[O,27r]'= 



,9i, . . . ,9j, . . . , 9k)+ 



+FN,k-i{di, ■ ■ ■ ,9i, ■ ■ ■ ,9j, . . . ,9k) 



- 9, 



9^$ 
92 $ 



I, ... ,9k) F/v,fc-i(6'i, . . . ,9i, . . . ,9j, . . . ,9k) 



+ 



''1, • • • , ^fc) -^A 



9i, ... ,9 j, ... ,9k 



{N - k)e^a^ ^ 
4 .^-^ 



^^2 V"i' • • • ' "K/ Af.fc-ll.C'i, . . . , i/j, . . . , 

— $(6'i, . . . ,9k) FN,k{9i, . . . ,9k) 



- ^0 



(i^x . . . d9k 
(27r)^- 



— U(9i,...A)6[0,27r]fe 



2 • • • ) ^fc) F^^kidi, . . . ,9k) — + 0{e] 



{2n) 



Now, hnking e and by fl4.12p . the 0(1) terms as — )■ 00 are on the first, second, fifth 
and sixth hues of the right-hand side. We can also drop the 0{e) and 0{e^) remainders. 
All the terms which are dropped are estimated by Ce||$||3oo uniformly in A^ and time, 
thanks to the property that F^ a; is a probability and to the assumptions of the theorem. 
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Therefore, F[k] satisfies 

d_ 
di 



J [0,2-^]^ 



d9i... d9k 



V 



E 

.i<j<k 



$(^1, . . . ,9k) ( ■ ■ ■ ,9i, . . . ,9j, . . . , 9k)+ 

d9i... d9k 



'(ei,...,efc)e[o,27r]fc 

+F[k-i]{9i, . . . ,9i, . . . ,9j, . . . ,9k) 5{9j — 9i) — 2F[k]{9i, . . .,9k 



" J(ei,...A)e[o,27r]fc d9i 



■ ■■ ,9k) F[k]i0i, ■ ■ ■ ,9k 



(27r)*^ 

d9i... d9k 

(27r)'= 



(8.13) 



Now, using Green's formula, Eq. f l8.13p appears as the weak form of f l4.13p . which ends 
the proof of Theorem 14.21 ■ 
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